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Abstract — This paper establishes methodology for the robust 
and stable design of infinite impulse response (IIR) digital 
filters using hybrid differential evolution method. Differential 
Evolution (DE) is undertaken as a global search technique 
and exploratory search is exploited as a local search technique. 
DE is a population based stochastic real parameter 
optimization technique relating to evolutionary computation, 
whose simple yet powerful and straight forward features make 
it very attractive for numerical optimization. Exploratory 
search aims to fine tune the solution locally in promising 
search area. This proposed DE method augments the capability 
to explore and exploit the search space locally as well globally 
to achieve the optimal filter design parameters by applying 
the opposition learning strategy and random migration. A 
multivariable optimization is employed as the design criterion 
to obtain the optimal stable IIR filter that minimizes the 
magnitude approximation error and ripple magnitude. DE 
method is implemented to design low-pass, high-pass, band- 
pass, and band-stop digital HR filters. The achieved design of 
IIR digital filters by applying DE method authenticates that 
its results are comparable to other algorithms and can be 
effectively applied for higher filter design. 

Index Terms — Digital IIR filters, Differential Evolution, 
Exploratory search algorithm, Multi-parameter, optimization, 
Opposition based learning, Random migration. 

I. Introduction 

The need for numerical optimization algorithms arises from 
almost every field of engineering, science, and business. This 
is because; an analytical optimal solution is difficult to ob- 
tain even for relatively simple application problems. A nu- 
merical algorithm is expected to perform the task of global 
optimization of a multimodal and non linear objective func- 
tions. However, an objective function possesses numerous 
local optima, which could trap the numerical algorithms. The 
possibility of failing to locate the desired global solution 
increases with the Increase in the problem dimension. The 
main challenge in designing infinite impulse response (IIR) 
digital filters is its nonlinearity and stability aspects in com- 
parison to finite impulse response (FIR) digital filter design. 
The digital filters can be implemented through either hard- 
ware or software. These are capable of processing both real- 
time and on-line signals. Image processing, speech synthe- 
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sis, secure communication, radar processing and biomedical 
etc are some of the areas where digital filters are useful. 

The design of digital IIR filters is achieved by either 
transformation or optimization technique. In the 
transformation technique, analog IIR filter is designed and 
then it is converted to digital IIR filter. By implementing 
transformation techniques [10], Butterworth, Chebyshev, and 
Elliptic functions, have been designed. Optimization methods 
have been applied whereby magnitude error, and ripple 
magnitudes of pass-band and stop-band are undertaken to 
measure performance of digital IIR filter design. Jiang et.al 
[30] has discussed the design of IIR digital filter by comprising 
stability constraint and has applied an iterative second-order 
cone programming method. Lightener et.al.[3] has discussed 
the simultaneous design in magnitude and group delay. A 
stability constraint with a prescribed pole radius has been 
derived from the argument principle of complex analysis for 
designing IIR digital filter. The semi-definite programming 
relaxation technique [31] has been implemented for design of 
IIR filter in a convex form. 

The design procedure is sequential and finds a feasible 
solution within a set of relaxed constraints. Normally 
conventional gradient-based design may easily get stuck in 
the local minima of error surface of IIR filters, due to its non- 
linear and multimodal nature. To overcome the shortcomings 
of gradient methods, researchers have applied nature based 
optimization algorithms such as genetic algorithms 
[4,6,8,9,1 1,15,16], particle swarm optimization (PSO) [17], 
seeker- optimization- algorithm -based evolutionary method 
[20], simulated annealing (SA) [12], Tabu search [19], ant 
colony optimization [18], immune algorithm [22] etc for the 
design of digital filters. 

Natural selection and genetics have been the basis for 
evolutionary algorithms and these are used as promising 
optimization techniques. Genetic algorithm (GA) also falls in 
the category of evolutionary algorithms. GA is capable of 
searching multidimensional and multimodal decision spaces. 
These are also useful to optimize complex and discontinuous 
functions [8] . The digital IIR filter can be structured as cascade, 
parallel, or lattice. Tang et al. [8] have applied genetic 
algorithms for designing the digital IIR filters. Genetic 
algorithms sometimes encounter very slow convergence. For 
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high dimensional optimization problems GA may trap in the 
local optima of objective function due to the existence of 
numerous local optima [7]. Tsai et al. [21] have put on the 
hybrid Taguchi genetic algorithm for design of optimal IIR 
digital filters. In hybrid Taguchi genetic algorithm, Taugchi 
has applied the combination of the traditional genetic 
algorithms .Consequently it becomes necessary to explore 
an algorithm that may provide global solution for the design 
of the optimal digital IIR filters. 

Tsai et al. [22] have suggested the integrating of the 
immune algorithm and the Taguchi method and named it as 
Taguchi-immune algorithm (TIA), for the design of digital IIR 
filter. Yu et.al. [24] have proposed cooperative co-evolutionary 
GA. During design of digital IIR filter, the lowest filter order 
has been found by considering the magnitude and the phase 
response. Here the digital IIR filter structure and its 
coefficients have been coded separately. The simulated 
annealing has been applied for the coefficient species so 
that diversity is maintained, for arriving at global minima [12], 
and it may require too many function evaluations. The seeker- 
optimization-algorithm is good at local convergence and it 
can be implemented for global minima [20]. It might often 
require too many objective function evaluations. The 
literature survey reveals that, there are various methods with 
which the optimization problems under different conditions 
can be addressed Optimization methods are classified based 
on the nature of objective function and type of search space. 
The evaluation of optimum solutions could be 
computationally expensive in IIR filter design problems due 
to time consuming computer simulation. Hence, a method is 
required, which is able to provide optimum solution within a 
given time budget. 

Kennedy and Ebehart [5] developed a global search 
technique and named it as particle swarm optimization (PSO). 
It is based on simple concept, easy to implement, 
computationally fast and has robust search ability. In PSO 
the social evolvement knowledge is simulated by searching 
the optimum through population evolution which may include 
fittest solutions. Despite many advantages PSO has some 
shortcomings because the convergence behavior of PSO 
depends upon its parameters. By chance PSO parameters 
have been chosen incorrectly, this may result in divergent 
particle trajectories which results trapping into local minima 
[27]. Further high-dimensional optimization problems, may 
suffer from the premature convergence problem. Till date, 
many efforts have been made by the researchers to improve 
the performance of PSO, either through mathematical analysis 
or by improving the convergence property [14, 32] . 

In the past, to tackle numerically complex computational 
problems, population based algorithms, like evolutionary 
algorithms (EAs), particle swarm optimization, and differential 
evolution, have been developed by the researchers. However, 
the performance of these algorithms may vary extensively 
from one problem to another and there is no single best 
algorithm for all problems. Practically, there is an inherent 
risk associated while choosing an algorithm for the given 
problem, as in prior no researcher knows which algorithm will 
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give optimal result for the problem. Further finding the 
optimum solution may not be practical or even possible for 
many real-world applications. The best solution obtained by 
a given iteration is often considered to be very important. In 
other words, an algorithm that performs extremely well on 
some problems but very poorly on other problems. So there 
is a more desirable that works well on a large variety of 
problems. 

These days' evolutionary algorithms comprising of GA, 
EP, ES, (DE) and so on are mostly used for solving complex 
real life computational problems. Differential Evolution differs 
from conventional GA in its use of perturbing vectors, which 
are the difference between two randomly chosen parameter 
vectors, a concept borrowed from the operators of Nelder 
and Mead's simplex optimization technique. DE algorithm 
was first introduced by Storn and Price in 1995 and was 
successfully applied in the optimization of some well-known 
nonlinear, non-differentiable, and non-convex functions. Das 
and Suganthan [33] provided an overview of the engineering 
applications that have benefited from the powerful nature of 
DE. Differential evolution algorithm has a number of 
significant advantages. DE has the ability to find the true 
global minimum regardless of the initial parameter values. It 
is fast and simple with regard to application. It requires few 
control parameters. It has parallel processing nature and fast 
convergence. DE is capable of providing multiple solutions 
in a single run. The method is effective on integer, discrete 
and mixed parameter optimization. DE has ability to find the 
optimal solution for a nonlinear constrained optimization 
problem with penalty functions. Although DE has many 
advantages explained above, this has a number of 
disadvantages. In DE, there exist many trial vectors 
generation strategies out of which a few may be suitable for 
solving a particular problem. Ten mutation strategies are 
undertaken to add a vector differential to population vector 
of individuals. Moreover, three crucial control parameters 
involved in differential evolution algorithm, i.e., population 
size, scaling factor, and crossover rate, may significantly 
influence the optimization performance [26]. Hybrid 
algorithms have the scope to overcome the shortcomings of 
the evolutionary techniques. 

The objective of this paper is to propose a hybrid DE 
method for the design of IIR digital filters that randomly 
explores the search space globally as well locally. The 
opposition employees for population initialization and for 
generation jumping The migration has been applied to 
maintain the diversity and search space exploration. Locally 
the solutions are fine tuned with exploratory search. The 
scaling factor is varied over generations by dynamic systems 
evidencing chaotic behavior. The values of the filter 
coefficients are optimized with the differential evolution to 
achieve magnitude error and ripple magnitude of pass-band 
and stop-band as objective functions for optimization 
problem. Here constraints are taken care of by applying 
exterior penalty method. 

The paper is organized in four sections. Section 2 
describes the IIR filter design problem statement. The 
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underlying mechanism and details regarding the differential 
evolutionary algorithm for designing the optimal digital IIR 
filters is described in Section 3. The performance of the 
proposed differential evaluation method has been evaluated 
and achieved results are compared with the design results 
by Tang etal. [8], Tsai etal. [21] andTsai and Horng [22] for 
the low pass, high pass, band pass, and band stop filters in 
section 4. Finally, in Section 5 the conclusions and discussions 
are outlined. 

U. IIR Filter Design Problem 

In the digital filter design, a set of filter coefficients is 
determined, which should meet the performance specifications 
of pass-band width and its corresponding gain, stop-band 
width and its attenuation, band edge frequencies, and 
tolerable peak ripple in the pass-band and stop-band. The 
realization of IIR digital filter is stated by the recursive 
equation [34]: 



N M 

y(n)=Yu x k u(n - k) -Yu x N+jy (n ~ J> 
k=0 j=\ 



(1) 



H(z)-- 



N 

*=0 



M 



(2) 



The design of digital filter design problem involves evaluation 
of a set of filter coefficients, Xk and xn+j which meet the 
performance indices. Several first-order and second-order 
sections are cascaded together [7, 8] for realizing IIR filters. 
The cascaded transfer function of IIR filter is denoted 

by H(co,x) , involving the filter coefficients like, poles and 
zeros. For all types of filter, the structure of cascading type 
digital IIR filter is stated as [4]. 

„, . (p, i+xye-f \ (» l+Xie-' a >+x i+ ie- 2 ' a > } 
where: 

l = 2N + 4(k-Y) + 2and vector X=[x, x, ...xj denotes the 

filter coefficients of dimension Sx 1 with S = 2N + 4M + 1 . 
In the IIR filter, the coefficients are optimized such that the 
approximation error function for magnitude is to be minimized. 
The magnitude response is specified at K equally spaced 
discrete frequency points in pass-band and stop-band. The 

absolute error is denoted by e(x) and is stated below: 



<?(*) = £ \H d (m ) - \H(co, ,x)\\ 



(4) 



Desired magnitude response, H d («, ) of IIR filter is given 
as: 



H d (co,)-- 



1, ;for (di e passband 



10, ;for cOjEstopband ^ 

The ripple magnitudes of pass-band and stop-band are given 
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by <5!(x)and 8 2 (x), respectively [2]. Ripple magnitudes for 
pass-band is given as: 

5l(A:)=max|//(«i,A:)|}-min|//(ffl 1 ',jr)|} passband 

coi coi ' {b) 

Ripple magnitudes for stop-band is stated as 

S 2 (x) = max|/7 e stopband 



coi 



(7) 



Stability constraints are included in the design of casual 
recursive filters, which are obtained by using the Jury method 
[1]. The multivariable constrained optimization problem is 
stated as below: 

Minimize f(x) = e(x) (8 a) 

Subject to the stability constraints:- 

l+x 2 , +1 >0 (i=l,2,....,A0 (8b) 

l-x 2l+1 >0 (j=l,2,....,A0 (8c) 

l-x, +3 >0 (l = 2N+4(k-l)+2,k=l,2,....,M) (8d) 

l+x, +2 +x, +3 >0(/=2Af+4(A:-l)+2,A:=l,2,....^f,) (8e) 

l-x, +2 +x, +3 >0(Z=2Af+4(A:-l)+2,A:=l,2,....,AO (8f) 

Scalar constrained optimization problem is converted into 
unconstrained multivariable optimization problem using 
penalty method. Augmented function is defined as 

A(x) = e(x) + r{P term ) (9) 

Where 

H 2 ^ 2 & 2 



M 



M 



+ ^(l+Jc;+2+Jc;+3) +^(l-^+2 + J?+3) 
k=\ k=l 



(10) 



r is a penalty parameter having large value. 

Bracket function for constraint given in (8b) is stated 
below: - 

l + x 2i+ i ;if(l + x 2i+1 \<0 



C + ^^lb"^' ';(fIl + "^>b(H) (ID 

Bracket function for constraint given by (8e) is stated 
below: - 

/ v_ Jl + x, +2 + x, +3 if (l + x, +2 + x, +3 )<0 

(i + x +2+ x )+3 ;- Q + (12) 



if (l + x, +2 + x, +3 ): 

Similarly bracket functions for other constraints given by 
(8c), (8d) and (8f) are undertaken. 

III. Differential Evolution for IIR Filter Design 

Differential Evolution (DE) is a population-based 
stochastic method. It is applied to minimize performance index. 
Differential evolution combines simple arithmetical operators 
with the classical operators of the recombination, mutation, 
and selection to evolve from a randomly generated starting 
population to a final solution [28]. The different variants of 
DE are classified using: DE/a/p/8. a indicates the method for 
selecting the parent chromosome that will form the base of 
the mutated vector. P indicates the number of difference 
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vectors used to perturb the base chromosome. 5 indicates 
the recombination mechanism used to create the offspring 
population. The bin acronym indicates that the recombination 
is controlled by a series of independent binomial experiments 
[33]. 

A Parameter Setup 

The user selects the key parameters that control the DE, 
i.e. population size (L), boundary constraints of optimization 

variables (S), mutation factor (f m ) , crossover rate (CR), and 
the stopping criterion of maximum number of iterations 

(generations) T max . X'j is the j"' element of S set of filter 

coefficients giving i"' individual of a population. The set of 
real IIR digital filter co-efficient (X) of all generators is 
represented as the population. For a system with S filter 
coefficients, the population is represented as a vector of 
length , S. If there are L members in the population, the 
complete population is represented as a matrix shown: 

B Initialization of an Individual Population 

Set generation t = 0. Initialize a population xjj (j=l,2, 

S; i = 1, 2, L) individuals with random values generated 
according to a uniform probability distribution in the S- 
dimensional problem space. Initialize the entire solution vector 
population within the given upper and lower limits of the 
search space. 

4 xf +rand()(xj ax -xf n ) 

(j=l,2,...,S; i=l,2,...,L) (13) 

where randf) is uniform random number between and 1 

The vector population may violate inequality constraints. 
This violation is corrected by fixing them either at lower or at 
upper limit. 

C. Opposition-Based Learning: 

Evolutionary optimization methods start with some initial 
random solutions and try to improve them toward some 
optimal solution(s). The computation time, among others, is 
related to the distance of these initial guesses from the optimal 
solution. It can improve the chance of starting with a better 
solution by simultaneously checking the opposite solution 
[23]. By doing this, the better one either guess or opposite 
guess can be chosen as an initial solution. As per the 
probability theory, 50% of the time, a guess is farther from 
the solution than its opposite guess [25]. Therefore, starting 
with the closer of the two guesses as judged by its objective 
function has the potential to accelerate convergence. The 
same approach can be applied not only to initial solutions 
but also continuously to each solution in the current 
population [25]. 

t _ mill max _ t 
X i+L,j ~ x j + x j x ij 

(j=l,2,...,S; i=l,2,...,L) 

(14) 
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where x™ n and ^™ ax are lower and upper limits of filter 
coefficients. 

D Evaluation of the Individual Population 

The goal is to minimize the objective function. The 

elements of parent/offspring Xy may violate constraint. A 

penalty term is introduced in the objective function to penalize 
its objective function value. Objective function is changed 
to the following generalized form: 

A i (x ij ) = e i (x ij ) + r(P term ) 

(j=l,2,..., S;i=l,2,,L) (15) 
where the penalty factor is given by (1 1) and (12). 

E Mutation Operation (Differential Operation) 

Mutation is an operation that adds a vector differential to 
a population vector of individuals. The mutation operation 
using the difference between two randomly selected 
individuals may cause the mutant individual to escape from 
the search domain So, ten variations mutation strategies [29] 



defined below are considered for study: 

Z ijl = P RJ + fm { x R 2 j ~ x R 3 j ) ( 1 6) 

Z ij2 = x Bj + fm [ X RJ ~ x R 2 j ) (17) 

Z \jh = x ij + fs [ x Bj ~ x ij )+ fm { X RJ ~ x R 2 j ) (18) 

Z ljA = x Bj + fm { X RJ + x R 2 j ~ x R 3 j ~ X RJ ) (19) 

Z lj5 = x R 5 j + fm { X RJ + x R 2 j ~ x R 3 j ~ X RJ ) (20) 

z U= x Bj+f m ( x Bj-4) (2D 

Z ljl = X Bj + fm { x Bj ~ x ij ~ X RJ ~ x R 2 j ) ( 22 ) 

Z ijS = x Bj + /b [ x Bj ~ P xij )+ fm { X RJ ~ x R 2 j ) (23) 

Z ij9 = x Bj + fm { x Bj + 4 ~ X RJ ~ x R 2 j ) (24) 

Z yiQ= x Bj + f m h t Bj~ x B~/) (25) 

(j=l,2,...,S;i=l,2,...,L) 



where t is the time (generation); R p R 2 and R 3 are mutually 
different integers t from the running index, i, randomly selected 
with uniform distribution. f,„(t)>0 is the mutation factor and 
is taken from the range [0.4, 1] using chaotic sequence. 
Better solution out of the ten mutation strategies based on 
minimum augmented objective function is selected. 

MZfj) = mm[M(Z! n ),A2(Zf j2 ), .... Am(Zf jW )\ (26) 

zj = [z/j, z? 2 , z' is J stands for the position of the ith individual 
of a mutant vector. 

F Recombination Operation 

Recombination is employed to generate a trial vector by 
replacing certain parameters of the target vector by the 
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corresponding parameters of a randomly generated donor 
vector. 

For each vector, Z ; t+1 , an index R 5 (i) e { 1, 2, S} is randomly 
chosen using a uniform distribution, and a trial 

vector,^' =\j%\u\i\...,U <?] 

jjM = fa # < CR) or (j = R 5 (ff) 

U * \P$ if(R 4 (j)>CR)or(j*R 5 (i)) 

(j=l,2, ...,S;i=l,2, ...,L) (27) 
where 

R 4 (j) is the jth evaluation of a uniform random number 
generation with [0, 1 ] . CR is the crossover or recombination 
rate in the range [0, 1]. 

Usually, the performance of a DE algorithm depends on three 
variables: the population size, the mutation factor f m (t) 
and the CR. 

G Selection Operation 

Selection is the procedure whereby better offspring are 

produced. To decide whether the vector U ' +1 should be a 
member of the population comprising the next generation, it 
is compared with the corresponding vector X' . Thus, if A 
denotes the objective function under minimization, then 

x t+\ = u = u - S) ' Jf A(t/ ' +1) < A{x 'i } 

V )x'jj (j=X2,..,S) otherwise 

(i= 1,2, ...,L) (28) 
In this case, the objective A of each trial vector Uff 1 is 

compared with that of its parent target vector 4j • If the 

augmented objective function, A of the target vector xfj is 

lower than that of the trial vector, the target is allowed to 
advance to the next generation. Otherwise, a trial vector 
replaces the target vector in the next generation. 

H. Exploratory Move 

In the exploratory move, the current solution denoting 
filter coefficients is perturbed in positive and negative 
directions along each variable one at a time and the best 
point is recorded. The current point is updated to the improved 
solution at the end of each filter coefficient perturbation. If 
the point found at the end of all filter coefficient perturbations 
is different from the original point, the exploratory move is a 
success; otherwise, the exploratory move is a failure. The 
filter coefficient x, is perturbed as follows: 



where 



xf=xf±S iU { (i=l2,...,S;j=l2,...,S) 



1 if i = j 
if i* j 



(29) 



(30) 



S denotes number of variables. 

The objective function denoted by A(xf)is calculated as 
follows [35] 



x° + A,. u t ; A(x° + A,. ) < A{x° ) 



-A,.«,. ;A(x 1 -A 1 . M/ )<A(x I ) 



;otherwise 



(31) 



where ( i = 1, 2,..., S ) and A; is random for global search and 

fixed for local search. 

The process is repeated till all the filter coefficients are 
explored and overall minimum is selected as new starting 
point for next iteration. 

J Random Migration Operator 

The population diversity and its exploration of the search 
space are rapidly decreased, and the clustered individuals 
cannot reproduce newly better individuals by mutation and 
crossover. In order to increase the exploration of the search 
space and decrease the selection pressure for a small 
population, it is randomly selected 0.2L individuals to start 
migration operation. The j' h gene of the i' h individual is 
randomly regenerated as follows [13]: 



x'f l + Ri(xf n -x'f l ) if8<-9- 



xlf+Riixf^-xft 1 ) if otherwise 



(32) 



where x' h J l is the best individual. R and 5 are uniform random 
number. 

K Verification of the Stopping Criterion 

Generation number is updated, ( t = t + 1). Procedure is 
repeated until a stopping criterion is met, usually a maximum 
number of iterations (generations), T max . The stopping 
criterion depends on the type of problem. 
L. ALGORITHM: IIR Filter Design Using differential 
Evolution Algorithm 

The search procedure of the proposed differential evolution 
method has been outlined below. 

1 Read data viz. Population size (L), mutation factor (f m ) , 
crossover rate (CR), and the stopping criterion of maximum 
number of iterations (generations) T max ., seed number, S e , 

xf n and x™ ax (i=l,2,..,S)etc. 

2. Generate an array of (SxL) size of uniform random numbers. 
FORj=l toL 
FORi=l toS 

3. 1 Generate the initial population individual, x~ using (13) 

3.2 Compute augmented objective function A t (xy) using 
(15). 

-,cACEEE 



©2013 ACEEE 
DOL01.IJSIR4.1.15 



Full Paper 



ACEEE Int. J. on Signal & Image Processing, Vol. 4, No. 1, Jan 2013 



3.2 Generate the initial population individual using opposition, 
xf +Lj using (14) 

3.3 Compute augmented objective function 

A i +L ( x i+Lj) using (14). 
END FOR 
END FOR 

4. Arrange Aj in ascending order and select first L population 

members out of 2L members. 

5. Set iteration counter, t = 

6. Increment the iteration counter, t = t + 1 

7. Select best member pbest and corresponding x' B j 
FORj = 1 to L 

8. 1 Generate an array of uniform random numbers and generate 
five different integer random numbers within the range of 1 
toL. 

8.2. Apply mutation operation to compute 
Z' ijk (i = 1,2,..., S;k = 1,2,.. .,10) using (16 to 25) 

8.3. Compute augmented objective function A ik (Zy k ) using 
(15). 

8.4. Find the mutant vector based on minimum augmented 
objective function using (26) to get A^Zy) and 

corresponding Z'j 

END FOR 
FORj=l toL 

9. 1 Generate arrays R4 and R5 of random numbers of size, 2L 

9.2. Apply recombination operation to compute 

C/fusing (27). 

9.3. Apply selection operation to compute variable 
x\J x (i = 1,2,..., S) using (28) and augmented objective 

function A t (x^ ) using (15). 

9.4. Apply exploratory move to improve the population 
by implementing algorithm 1 . 

9.5. Apply random migration to compute variable using 
(32) and augmented objective function using (15). 

END FOR 

10. IF(t<T ) THEN GOTO 6 

v max' 

11. STOP. 

IV. Design of IIR Filters And Comparisons 

The design of cascaded digital IIR filter has been 
implemented by evaluating filter coefficients using differential 
evolution. Low-pass (LP), high-pass (HP), band-pass (BP) 
and band-stop (BS) filters have been considered for the 
design. To design digital IIR filter, 200 equally spaced points 
are set within the frequency domain [cu] , 



A. Low-Pass (LP) Filter: 

The predefined range of pass-band and stop-band are 
taken as 0<ffl<0.2;r and 0.3;r <a< n , respectively. The 
value of N and M has been taken as 1 . The maximum number 
of iterations has been taken as 200 for DE. The maximum 
migration value is 50. The rate of opposition varies between 
and 1 and has been taken as 0.6. The mutation ratio f and 

m 

crossover ratio CR, has been taken as 0.85 and 0.25, 
respectively. Exploratory move is repeated 20 times. 
The low pass IIR filter model designed by the DE approach 
is given below in equation (33). 



H LP (z) = 0.036514 x 



(z + 0.8925 1 )(z 2 -0.27986 z + 0.948 10 ) 
(z -0.66080 )(z 2 -1 .39399 z +0.72953 ) 



(33) 




Figure 1: Magnitude with respect to band-pass ripples and band-stop 
ripples graph of Low-Pass Filter 

The variation in magnitude with respect to the ripples of 
pass-band and stop-band for LP filter has been drawn from 
100 iterations in Fig. 1 . The minimum value of magnitude is 
stabilized after 80 generations as shown in Fig. 2. Out of 100 
iterations run, the obtained best values for LP filter design 
are given in Table I. The pass-band and stop-band 
performance has been obtained from (6) and (7) respectively 
and are given in Table I. 



r.i* 




Figure 2: Magnitude with respect to Iteration Graph of LP Filter 

The magnitude response of LP filter with respect to nor- 
malized frequency has been shown in Fig. 3. The pole-zero 
graph of low-pass filter has been depicted in Fig. 4 which 
shows the stability of the designed filter. The magnitude sta- 
bilizes after 80 generations. 
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Table I. Design Result For LP Filter 



Method 


Magnitude 
Error 


Pa35-band performance 


Stop-band 
performance 


DE 


3i014 


C.3333e<|Jf(e^))<l.<H9 
(0.1353) 


|ff(s.*)|<O.l505 
(0.1.505) 


HGA 
[3] 




C.SS70< |Sf(s^)|< 1.009 
(01139) 


|H(^;|£C.1SC2 
(0.1302) 


HTGA 
[21] 


42511 


C.SCCC4i|H<s, i ^<1000 
(0.0996) 


|H(s^)|<0.1247 
(0.1247) 


ha 

[22] 


3.S157 


O^M^IflV^ 1.000 
(0.1036) 


|ff(s/^))<0.1633 
(0.163S) 




0.4 0.6 
normalized frequency 



0.3 



1 



Figure 



3: Magnitude with respect to normal frequency response 
of LP filter 

r" 




u 

Seal li.i1 

Figure-4: Pole-Zero graph of low pass filter 

B. High Pass (HP) Filter: 

The predefined range of pass-band and stop-band are 
taken as 0M<o}<tt and 0<»<0.77r , respectively. The N and 
M value has been taken as 1. The maximum number of 
iterations has been taken as 200 for DE. The maximum 
migration value is 50. The rate of opposition varies between 
and 1 and has been taken as 0.6. The mutation ratio f and 
crossover ratio CR, has been taken as 0.85 and 0.25, 
respectively. Exploratory move is repeated 20 times. The high- 
©2013ACEEE 
DOI:01.IJSLP4.1.15 



pass IIR filter model designed by the DE approach is given 
below in equation (34). The obtained best values for HP filter 
design are given in Table II. The magnitude response of HP 
filter has been shown in Fig. 5 and the pole-zero graph of 
high-pass filter has been depicted in Fig. 6 which shows the 
stability of the designed filter. 



H HP (z)= 0.015683 > 
' z 2 + 0.56873 



z- 0.90131 ] 
z + 0.75157 J 
z + 0.79047 "| (34) 



z 2 + 1.63458 z + 0.78028 



Table II. Design Results For HP Filter 



Method 


Magnitude 
Error 


Pass -band 
perform, snee 


Stop-baud 
performance 


DE 


2.8960 


C.3S553<|H(s>)|<1.014 
(0.11889) 


|ff(s,^))<ftllO0 
(0.1100) 


HGA 
[8] 


14.5078 


C.9224<|ff(.s-*}|< 1.CC1 
(0.0779) 


ff(s.*)|<0.1B19 

(0.1819) 


HTGA 
[21] 


4.8372 


C.946C< Jf(s;^<l.C0C 
(0.0540) 


\H(s.:*p«.U57 
(0.1457) 


TIA 
[22] 


4.1819 


<s.m9<ft(e;«pi.m> 
(0.0771) 


|H(a : *;|<C.1424 

(0.1424) 




Figure 6: Pole-Zero Graph of High Pass Filter 
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C. Band Pass (BP) Filter: 

The predefined range of pass-band and stop-band are 
taken as QAn < co < 0.6k and 0<a<0.25n: , 0.15<m<n: , 
respectively. The value of N and M has been taken as and 
3. The maximum number of iterations has been taken as 200 
for DE. The maximum migration value is 50. The rate of 
opposition varies between and 1 and has been taken as 0.6. 
The mutation ratio f and crossover ratio CR, has been taken 
as 0.85 and 0.25, respectively. Exploratory move is repeated 
20 times. The band pass IIR filter model designed by the DE 
approach is given below in equation (35). 



H B p(z) = 0.019873 



-0.00063 



z z ■ 



-0.00260 z- 1.36653 



z 2 + 0.00476 z + 0.55156 
-1.19862 ^ 



z 2 - 0.62427 z + 0.77647 
z 2 +0.00074 z - .922127 
z 2 + 0.62495 z + 0.77346 



(35) 

The obtained best values for BP filter design are given in 
Table III. The magnitude response of BP filter has been shown 
in Fig. 7 and the pole-zero graph of band- pass filter has been 
depicted in Fig. 8 which shows the stability of the designed 
filter. 

Table III. Design Results For Bp Filter 



Method 


Magnitude 
Error 


Pas s -band 
perform auoe 


Stop-band 
perform auce 


DE 


1.25 SO 


(0.0226) 


|H<a/»]|<0.05n 
(0.0512) 


EGA 
[S] 


5.2165 


0SS;r^|H{e.*:|<1.000 
(0.1044) 


|ff{e;»)|<0.1772 
[0.1772) 


HTGA 
[21] 


1941S 


(0.0234) 


|H(af*3|< 0.0711 

(0.0711) 


TIA 

[22] 


1 5204 


0S6Sl<|ff<s/*)|< 1.000 

(0.0319) 


\H{b;°$< 0.0679 

(0.0679) 




OA 0.6 
Normalised frequency 

Figure 7: Magnitude with respect to normal frequency response of 
BP Filter 
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D. Band Stop (BS) Filter: 

The predefined range of pass-band and stop-band are 
taken aso<w<0.25ff 075<w<tt and 0.4^ <w<0.6^ , 
respectively. The N and M values has been taken as and 2. 
The maximum number of iterations has been taken as 200 for 
DE. The maximum migration value is 50. The rate of opposition 
varies between and 1 and has been taken as 0.6. The 
mutation ratio f and crossover ratio CR, has been taken as 

m 

0.85 and 0.25, respectively. Exploratory move is repeated 20 
times. The band-stop IIR filter model designed by the DE 
approach is given below in equation. (36) and the best results 
obtained are given in Table IV. 




■05 0.5 

F»a par 

Figure 8: Pole-zero graph of Band-Pass 



H BS (z) = 0.415583 



z 2 + 0.12044 z + 0.99614 
z 2 + 0.58857 z + 0.50540 
0.70514 z + 0.98496 



-0.95940 z + 0.54561 



Table IV. Design Results For Bs Filter 



(36) 



Method 


Magnitude 
Error 


Pass-band 
perform moe 


Stop-band 
perfoimmoe 


DE 


3.13S5 


0.92 62 Zp(s;-°$< 1.010 

(o.om) 


\H(^;^< 0.1632 
(0.1632) 


HGA [S] 


6.6072 


C.3WC<|ff(s^)|<1.000 
(O.IOSO) 


|H(s..-^|<0.1726 
(0.1726) 


HTGA 
[21] 


45504 


093 mp(s;"p 1.000 
(0.0437) 


p&°>)\< 0.1013 
(01013) 


TIA [22] 


14750 


0.92 5 9<p(s^;\< 1.000 
(00741) 


\H(s;K\< 01273 
(0.1278) 



The magnitude response of BS filter has been shown in 
Fig. 9 and the pole-zero graph of band-stop filter has been 
depicted in Fig. 10 which shows the stability of the designed 
filter. 
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The Numerical results are compared with the filter design 
results given by Tang et al. [8], HTGA [21] and TIA [22] for 
the LP, HP, BP, and BS filters. To design the digital IIR filters, 
the augmented objective function with the stability con- 
straints, given by (9) is minimized. 1 

That is for the digital IIR filters designed under the pro- 
posed method gives better performances than the GA-based 
approach given in HGA [8] , HTGA [21] and TIA [22] ap- 
proaches. In case of Low-pass filter, it has been observed 
from Fig. 1 that magnitude of ripples in band-stop increases if 
the magnitude of ripples in band-pass decreases and vice 
versa. 




T) 0.2 Q.4 OS 1 

Nortnaloed rrequenty 

Figure 9: Magnitude with respect to normal frequency response of 
BS Filter 
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Figure 10: Pole-zero plot of band-stop filter 



Table V. Maximum, Minimum, Average Values And Standard Deviation of 
Magnitude Error 



Fdter 
Type 


Maximum 
Magnitude 
Error 


Minimum 
Magnitude 
Error 


Average 
Mapitude 
Error 


Standard 
Deviation of 
Magnitude 
Error 


Lew 
Pass 


6211293 


3.434767 


4.54022 


0.669059 


High 
Pass 


6.234363 


2.7S26S6 


4.31694 


0.997969 


Band 
Stop 


5.261713 


3.073425 


4.039234 


0.497461 


Band 
Pass 


1.639933 


1 2343714 


1.42674 


0.I1S4S3 
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Similar trend has been observed for HP, BP and BS filters. 
Moreover, 100 independent trial runs have been given with 
random seed numbers to have different uniform random 
numbers for studying the variation in the magnitude. 

The maximum value of magnitude error, minimum value of 
magnitude error, average value of magnitude error and 
standard deviation in magnitude error are given in Table V 
The results obtained depict that the standard deviation (SD) 
is less than 1.0 which shows the robustness of the designed 
system. 

Conclusions 

This paper proposes hybrid differential evolution method 
for the design of digital IIR filters whereby locally fine tuned 
by exploratory search method. As shown through simulation 
results, the DE method works well with an arbitrary random 
initialization and it satisfies prescribed amplitude 
specifications consistently. Therefore, the proposed 
algorithm is a useful tool for the design of IIR filters. 

On the basis of above results obtained for the design of 
digital IIR filter, it can be concluded that the proposed DE 
method hybridized with exploratory search is superior to the 
GA-based methods. Further, the proposed DE approach for 
the design of digital IIR filers allows each filter, whether it is 
LP, HP, BP, or BS filter, to be independently designed. The 
proposed DE is very much feasible to design the digital IIR 
filters, particularly with the complicated constraints. 
Parameters tuning still is the potential area for further research. 
The unique combination of exploration search and global 
search optimization method that is DE provided by the two 
types of algorithms yields a powerful option for the design 
of IIR filters. 

References 

[1] I. Jury, Theory and Application of the Z-Transform Method, 
New York: Wiley, 1964. 

[2] M.R. Lightner and S.W. Director, Multiple criterion 
optimization for the design of electronic circuits, IEEE 
Transactions on Circuits and Systems, vol. 28, no. 3,pp. 169- 
179, March 1981 

[3] G. Cortelazzo and M.R. Lightner, Simultaneous design in both 
magnitude and group-delay of IIR and FIR filters based on 
multiple criterion optimization, IEEE Transactions on 
Acoustics, Speech, and Signal Processing, vol. 32, no. 5 pp. 
949-967, 1984. 

[4] S .C . Ng, C . Y. Chung, S .H . Leung, and A. Luk, Fast convergent 
genetic search for adaptive LIR filtering, Proceeding IEEE 
International Conference on Acoustic, Speech and Signal 
Processing, Adelaide, pp. 105-108, 1994 

[5] J. Kennedy and R. Ebehart, Particle swarm optimization 
Proceedings of LEEE International Conference Neural Network, 
Perth Australia, No 4, pp.1942-1948, 1995 

[6] J.H. Li and EL. Yin, Genetic optimization algorithm for designing 
IIR digital filters, Journal of China Institute of 
Communications, vol. 17, pp. 1-7, 1996. 

[7] J.M. Renders and S.P. Flasse, Hybrid methods using genetic 
algorithms for global optimization, IEEE Transactions on 
Systems, Man, and Cybernetics-Part B: vol. 26, no. 2, pp. 

vr ACEEE 



Full Paper 



ACEEE Int. J. on Signal & Image Processing, Vol. 4, No. 1, Jan 2013 



243-258, April 1996. 

[8] K.S. Tang, K.F. Man, S. Kwong and Z.F. Liu, Design and 
optimization of IIR filter structure using hierarchical genetic 
algorithms, IEEE Transactions on Industrial Electronics, vol. 
45, no. 3, pp. 481-487, June 1998. 

[9] S.P. Harris and E.C. Ifeachor, Automatic design of frequency 
sampling filters by hybrid genetic algorithm techniques, IEEE 
Transactions on Signal Processing, vol. 46, no. 12, pp. 3304- 
3314, December 1998. 

[10] AV. Oppenheim, R.W. Schafer and J.R. Buck, -Discrete-Time- 
Signal-Processing, Prentice Hall, NJ, Englewood Cliffs, 1999. 

[11] K. Uesaka and M. Kawamata, Synthesis of low-sensitivity 
second order digital filters using genetic programming with 
automatically defined functions, IEEE Signal Processing Letters, 
vol. 7, pp. 83-85, April 2000. 

[12] S. Chen, R. H. Istepanian, and B. L. Luk, Digital IIR filter 
design using adaptive simulated annealing, of Digital Signal 
Processing, vol. 11, no. 3, pp. 241-251, July 2001. 

[13] G. Vanuytsel, P. Boets, L.V. Biesen and S. Temmerman, Efficient 
hybrid optimization of fixed-point cascaded IIR filter 
coefficients, Proceedings: IEEE International Conference on 
Instrumentation and Measurement Technology, Anchorage, 
AK, pp.793-797, 2002 
[14] M. Clerc and J. Kennedy, The particle swarm-explosion, 
stability, and convergence in a multi-dimensional complex 
space, IEEE Transactions on Evolutionary Computation 2002; 
vol. 6, no.l,pp 58-73, 2002 

[15] G.X. Zhang, W.D. Jin, and F. Jin, Multi-criterion satisfactory 
optimization method for designing IIR digital filters, 
Proceedings: International Conference on Communication 
Technology, Beijing, China, pp. 1484-1490, 2003 

[16] Y. Liu, S.Q. Shao, H. Zhao, X.F Liao, and LB. Yu, An 
application of genetic algorithms with guiding strategy in the 
design of digital filters, Proceeding: International Conference 
on Communication, Circuits Systems, Chengdu, China, pp. 
1141-1145,2004 

[17] J. Sun, W.B. Xu, and B. Feng, A global search strategy of 
quantum behaved particle swarm optimization, Proceeding 
IEEE Conference Cybernetics and Intelligent Systems, 
Singapore, pp. 111-116, 2004 

[1 8] N. Karaboga, A. Kalinli, and D. Karaboga, Designing IIR filters 
using ant colony optimization algorithm, Journal of 
Engineering Applications of Artificial Intelligence, vol. 17, no. 
3, pp. 301-309, April 2004. 

[19] A. Kalinli and N. Karaboga, A new method for adaptive IIR 
filter design based on Tabu search algorithm, International 
Journal of Electronics and Communication (AEU), vol. 59, 
no. 2, pp. 111-117, 2005. 

[20] Chaohua Dai, Weirong Chen and Yunfang Zhu, Seeker 
optimization algorithm for digital IIR filter design, IEEE 
Transactions on Industrial Electronics, vol. 57, no. 5, pp. 
1710-1718, May 2006. 

[21] Jinn-TsongTsai, Jyh-Horng Chou and Tung-Kuan Liu, Optimal 
design of digital IIR filters by using Hybrid Taguchi Genetic 
Algorithm, IEEE Transactions on Industrial Electronics, vol. 
53, no. 3, pp. 867-879, June 2006. 



[22] Jinn-Tsong Tsai and Jyh-Horng, Optimal design of digital IIR 
filters by using an improved immune algorithm, IEEE 
Transactions on Signal Processing, vol. 54, no. 12, pp. 4582- 
4596, December 2006. 

[23] H.R.Tizhoosh, Opposition based reinforcement learning, J 
advanced computation, Intel intelligent Inform, vol 10, no. 3, 
pp 578-585, 2006.1 

[24] Y. Yu and Y. Xinjie, Cooperative co-evolutionary genetic 
algorithm for digital IIR filter design, IEEE Transactions on 
Industrial Electronics, vol. 54, no. 3, pp. 1311-1318, June 
2007. 

[25] Rahnamayan, H.R. Tizhoosh and M.A. Salama, Opposition 
based differential evolution, IEEE Transactions on 
evolutionary computations, vol 12, no.l, pp 64-79, Feb 2008 

[26] Nasimul Noman and Hifoshi Iba, accelerating differential 
evolution using an adaptive local search, IEEE Transactions 
on evolutionary computation, vol 12, no l,pp 107-125, Feb 
2008 

[27] M. Omran, A. Engelbrecht and A. Salman, Bare bones 
differential evolution, European Journal of Operational 
Research; vol. 196, no.l, pp.128-139, 2009 

[28] A.K. Qin, V.L. Huang and P.N. Sugathan, Differential evaluation 
algorithm with strategy adapter for global numerical 
optimization, IEEE Transactions on Evolutionary 
Computation, vol 13, no. 2,pp 398-417,April 2009 

[29] Sum-Im, GA. Taylor, M.R. Irving and YH. Song, Differential 
evolution algorithm for static multistage transmission 
expansion planning, IET generation, distribution, vol 3(4), pp 
365-384, 2009 

[30] Aimin Jiang and Hon Keung Kwan, IIR digital filter design 
with new stability constraint based on Argument Principle, 
IEEE Transactions on Circuit and Systems-I, vol. 56, no. 3, 
pp. 583-593, March 2010. 
[31] Aimin Jiang and Hon Keung Kwan, Minimax design of IIR 
digital filters using SDP relaxation technique, IEEE 
Transactions on Circuits and Systems-I, vol. 57, no. 2, February 
2010. 

[32] Jun Sun, Wei Fang and Wenbo Xu, AQuantum-Behaved Particle 
Swarm Optimization with diversity-guided mutation for the 
design of two-dimensional IIR digital filters, IEEE Transactions 
on Circuits and Systems-II, vol. 57, no. 2, pp. 141-145, 
February 2010. 

[33] S. Das and P.N. Suganthan, Differential evolution: A survey of 
the state-of-the-art, IEEE Transactions on Evolutionary 
Computation, vol. 15, no. 1, pp. 4-31, February, 2011 

[34] E.C. Ifeachor and B.W. Jervis, Digital signal processing, a 
practical approach, 2nd edition, Pearson Education, Singapore, 
2003 

[35] J.S. Dhillon 1 , J.S. Dhillon 2 and D.P. Kothari, Economic- 
emission dispatch using binary successive approximation based 
evolutionary search, IET Proceedings-Generation, 
Transmission and Distribution, vol.3, no.l, January 2009, pp. 
1-16 



©2013 ACEEE 
DOI:01.IJSIP.4.1.15 



10 



-v'c ACEEE 



